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L  Scientific  Progress  and  accomplishments: 

The  major  accomplishments  in  this  work  can  be  put  in  three  categories  : 

1.  A  computer-based  approach  for  checking  die  quality  of  die  error  estimators, 

2.  A  practical  method  for  reliable  estimation  of  the  error  in  the  quantity  of  interest, 

3.  An  innovative  approach  for  adaptive  mesh-refinement  to  obtain  the  desired  quantity  of  interest  with  prescribed  accuracy. 

The  theoretical  foundations  for  the  above  accomplishments  and  die  significance  of  these  methods  in  practical  computations  are 
given  in  [1-22].  The  above  accomplishments  provide  a  strong  basis  towards  developing  robust  algorithms  for  the  simulation  of  the 
crack  propagation  using  adaptive  hp  finite  element  analysis. 

In  [1-3]  we  have  shown  that  the  error  in  the  finite  element  solution  at  a  point  can  be  split  into  two  parts  -  the  local  error  and  the 
pollution  error .  In  [3-14]  we  have  shown  that  standard  error  indicators  which  are  employed  in  practice  estimate  only  the  local  error 
and  hence  they  will  severely  underestimate  the  error  in  the  finite  element  solution  in  certain  cases.  We  have  developed  a  new  method 
for  the  a-posteriori  estimation  of  the  pollution  error  by  using  die  standard  error  indicators  for  the  solution  of  an  auxiliary  problem. 
This  method  can  be  implemented  in  any  commercial  finite  element  program  which  has  the  capability  to  solve  multiple  loads  and 
provide  an  estimate  of  die  global  energy  norm  of  the  error.  Thus  we  have  been  able  to  estimate  the  error  in  the  stress  in  the  “hot¬ 
spots”  as  demonstrated  in  [3,12,14,15,21].  We  have  also  shown  in  [3]  that  in  order  to  obtain  reliable  estimates  of  the  pollution  error, 
it  is  essential  to  employ  error  indicators  which  are  robust  in  the  interior  of  the  mesh,  at  the  boundary  and  in  the  vicinity  of  the 
comers. 

In  [1,2,4,9,18]  we  developed  a  computer-based  approach  for  checking  the  quality  of  die  error  indicators  for  the  class  of  meshes 
and  materials  employed  in  practice.  The  quality  of  the  error  indicators  is  often  measured  by  ejfectivity  index ,  which  cannot  be 
computed  in  practice.  However,  it  is  possible  to  determine  the  range  of  the  variation  of  the  effectivity  index  for  a  class  of  meshes, 
materials  and  types  of  problems.  We  introduced  the  concept  of  robustness-index  which  is  a  measure  of  the  range  of  the  effectivity 
index  of  an  error  indicator  for  a  class  of  meshes,  materials  and  types  of  problems  and  we  have  given  a  practical  method  for 
computing  the  robustness-index  for  any  error  indicator  which  maybe  also  given  in  the  form  of  a  “black-box”  in  a  commercial  finite 
element  program.  Based  on  this  approach  we  could  check  the  quality  of  several  error  estimators  which  are  used  in  practice  and  we 
have  also  suggested  modifications  for  the  existing  estimators  to  obtain  better  performance.  We  concluded  that  the  estimator  based  on 
least-squares  fit  of  the  stresses  (also  called  ZZSPR  estimator)  is  the  most  robust  estimator  for  meshes  of  triangles  and 
quandrilaterals  with  straight  edges.  We  have  also  concluded  that  the  commonly  used  explicit  estimator  can  lead  severe 
overestimation  of  the  energy  norm  of  the  error  for  particular  types  of  meshes  of  triangles. 

In  [14,15,21]  we  gave  a  new  approach  for  the  a-posteriori  estimation  of  the  error  in  the  quantity  of  interest  and  we  also  gave 
algorithms  for  the  adaptive  mesh-refinement  to  obtain  die  stress-intensity  factors,  for  a  body  with  several  cracks  and  reentrant 
comers,  within  the  prescribed  tolerance.  The  a-posteriori  estimation  of  the  error  in  die  quantity  of  interest  is  based  on  the 
construction  of  an  auxiliary  solution  which  is  computed  by  using  an  extra  right-hand  side  in  the  global  solution.  We  gave  several 
examples  to  illustrate  that  die  adaptive  mesh-refinement  based  on  controlling  the  global  energy  norm  of  the  error  does  not  guarantee 
the  prescribed  tolerance  for  the  stress-intensity  factors.  The  alogorithms  developed  in  [6,14,15,21]  lead  to  meshes  with  fewer 
degrees  of  freedom  than  those  obtained  from  the  ones  based  on  controlling  the  global  energy  norm  of  the  error  and  moreover,  we 
can  guarantee  the  accuracy  of  the  stress-intensity  factors  for  all  the  required  cracks  in  the  domain. 
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